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Abstract 

We present a unified geometric framework for describing both the Lagrangian and Hamil- 
tonian formalisms of regular and non-regular time-dependent mechanical systems, which is 
based on the approach of Skinner and Rusk [18) . The dynamical equations of motion and 
their compatibility and consistency are carefully studied, making clear that all the character- 
istics of the Lagrangian and the Hamiltonian formalisms are recovered in this formulation. 
As an example, it is studied a semidiscretization of the nonlinear wave equation proving the 
applicability of the proposed formalism. 

Key words: Lagrangian and Hamiltonian formalisms; autonomous mechanics, symplectic 
and presymplectic manifolds. 
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1 Introduction 

In 1983 Skinner and Rusk introduced a representation of the dynamics of an autonomous me- 
chanical system which combines the Lagrangian and Hamiltonian features [18]. The aim of this 
formulation was to obtain a common framework for both regular and singular dynamics, ob- 
taining simultaneously the Hamiltonian and Lagrangian formulations of the dynamics. Over the 
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years, however, Skinner and Rusk's framework was extended in many directions. So, Cantrijn et 
al [2] extended this formalism for explicit time-dependent systems using a jet bundle language. 
In [6] an extension of this formalism to other kinds of more general time-dependent singular dif- 
ferential equations was given. Cortes et al [3J used the Skinner and Rusk formalism to consider 
vakonomic mechanics and the comparison between the solutions of vakonomic and nonholonomic 
mechanics. Finally, in [1JE1IT5] the Skinner-Rusk model was developed for classical field theories. 

The aim of this paper is to continue the study of the the Skinner-Rusk formalism for time 
dependent mechanical systems (Section [3]), now, carefully studying the dynamical equations of 
motion and the submanifolds where they are consistent, and showing how the Lagrangian and 
Hamiltonian descriptions are recovered from this unified framework (Sections 1415ft . 

The case of field theories was independently developed in [HE], and improves the construction 
given in [2] , as it is discussed in Section 

As a new application, we analyze the case of semidiscretizations of field theories in Sectional 
These methods are designed by numerical schemes that respect physical principles preserved by 
the continuous systems, specially those described by partial differential equations (PDEs). In this 
case, there are not only a time dependence (as in ordinary differential equations) but also posses 
an spatial dependence. Many integration methods, in particular in Hamiltonian dynamics, starts 
by discretizing the spatial structure (spatial truncation) obtaining a finite dimensional system of 
ordinary differential equations (ODEs) retaining some physical properties of the original system 
(see [8]). For simplicity, we restrict ourselves to a particular semidiscretization of the nonlinear 
wave equation \10\ [T3] obtaining a unique solution of the dynamics on the secondary constraint 
submanifold. 

All the manifolds are real, second countable and C°°. The maps are assumed to be C°° . Sum 
over repeated indices is understood. 

2 Non- autonomous Lagrangian and Hamiltonian systems 

(See [U El [HI HU [16] for more details). In the jet bundle description of non-autonomous dy- 
namical systems, the configuration bundle is tt: E >• R, where E is a (n + l)-dimensional 

differentiable manifold endowed with local coordinates (t, q l ), and R has t as a global coordinate. 
The jet bundle of local sections of it, J 1 -^, is the velocity phase space of the system, with natural 
coordinates (t,q l ,v % ), adapted to the bundle tt: E »- R, and natural projections are 

tt 1 : JV *~ E , tt 1 : J 1 tt »- R . 

(If E = R x Q, where Q is a n-dimensional differentiable manifold, then J 1 tt ~ R x TQ). 

A Lagrangian density C G ^(J 1 ^) is a 7r 1 -semibasic 1-form on J lr rr, and it is usually written 
as £ = Ldt, where L G C 00 ^ 1 ^) is the Lagrangian function determined by C. Throughout this 
paper we denote by dt the volume form in R, and its pull-backs to all the manifolds. 

The Poincare-Cartan forms associated with the Lagrangian density C are defined using the 
vertical endomorphism V of the bundle J 1 tt (see [5[ I17j) 

e c = i{v)dc + c e ; n c = -do c g n 2 (j\) . 

A Lagrangian C is regular if £lc has maximal rank; elsewhere C is singular. In natural coordinates 
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d d 

we have V = (do* — v l dt) ® — - ® —, and 

ov l at 

^ - w*t -{&->)« 

= - 7n ^ i dviAd q i - 7 — n d q iAdq i 
ov-'ov 1 oq J ov l 

+ _&L_ v i dv j Adt+ ( d 2 L vi _ dL_ + d 2 L \ d j a dt 
dvidv 1 \dqidv i dqi dtdvi J 

The regularity condition is equivalent to det ( A y) ) ^ 0, for every y 6 J lr K. Geometrically, 

\ov l ov3 J 

C is regular if and only if (Oc,dt) is a cosymplectic structure on J 1 ^. This means that Qc an d 
dt are closed and O,^ A dt is a volume form (see 

The Lagrangian problem consists in finding sections (j>: R »- E of ir, characterized by 

(j'V)* i(-X")n£ =0 , for every X G £( J 1 ^) 

where j 1 ^ : R >■ J 1 ^ is the 1-jet extension of <f>. In natural coordinates, if </>(t) = (t,0*(t)), 

this condition is equivalent to demanding that <j) satisfies the Euler-Lagrange equations 



8L_ 

dq 1 



d_ f dL 

iH ~ ~dt I cV 



= , (for i = 1, ...,n) 

j 1 ^ 



where j l 4>{t) = (*j ^(t), 4> l {t)). Assuming that these sections are integral curves of vector fields 
in J 1 ^ the corresponding equations for these vector fields are 

i(x c )n c = o , i(x c )dt = i (i) 

where Xq_ G X(J 1 ir) is holonomic (recall that a vector field in J 1 ^ is said to be holonomic, or also 
a second order differential equation (SODE for simplicity), if its integral curves are holonomic; 

that is, canonical liftings of sections (p: R E). In the regular case, there is a unique solution 

to these equations. In the singular case the existence of a solution is not assured, except perhaps 
on some submanifold (or subset) of J^ir, where the solution is not unique, in general. 

Consider now the extended momentum phase space T*E, and the restricted momentum 
phase space which is defined by J 1 -^* = T*S/7r*T*R. Local coordinates in these manifolds are 
{t,q l ,p,Pi) and (t,q l ,pi), respectively. Then, the following natural projections are 

t^JV ^E , r 1 = vr o t 1 : jV ^ R , /z: T*E jV , p:T*E *-R. 



Let 9 G Q l (T*E) and n = -d@ G Q 2 (T*E) be the canonical forms of T*E whose local 
expressions are 

6 = pidq 1 + pdt , = dq 1 A dpi + dt A dp . 

(In the particular case E = R x Q, we have T*E ~ R x R* x T*Q, and JV ~ R x T*Q 

and introducing the projections prf. T*(R x Q) R x R*, pr 2 : T*(R x Q) *- T*Q, we 

have 6 = prl@ m + pr^&Q and = pr*0 R + pr|OQ; where K = -d6 K G i? 2 (R x R*) and 
Q.q = -dQ Q G Q 2 {T*Q) denote the natural symplectic forms of R x R* and T*Q). 

Being 0£ G f2 1 (J 1 7r) 7T -semibasic, we have a natural map T X: J 1 tt >■ T*E, given by 

FC{y) = &dy) (2) 
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which is called the extended Legendre map associated to the Lagrangian density C The re- 
stricted Legendre map is TL = [i o T X: J l n *- J 1 it*. Their local expressions are 

- — ■ — ■* * . . - — ■ — ■* (jTj - — ■ — ■* ■ oL 

TLt = t , TLq % = q l , TL Pi = — , TCp = L-v t — i 

o T 

TCt = t , Tl?j = q i , TC* Pi =— i 

or, in othj* words, TC{t, ?W) = (t, q i ,L - v*§k, f£) and F£(t, q\ q l ) = (t, q\ g). Moreover, 
we have TC = 0£, and TC VL = Sic- 

The Lagrangian C is regular if, and only if, J-C is a local diffeomorphism. As a particular 
case, £ is a hyper-regular Lagrangian if J-C is a global diffeomorphism. 

If C is a hyper-regular Lagrangian, then V = TC{J 1 tt) is a 1-codimensional, /i-transverse 
embedded submanifold of T*^, with natural embedding jo: V ^ T*E, which is diffeomorphic 
to J 1 -^*. This diffeomorphism is the inverse of p, restricted to V, and also coincides with the 
map h = J-CoJ r C~ 1 , when it is restricted onto its image (which is just V). This map h is called 
a Hamiltonian section, and is used to construct the Hamilton-Cartan forms in J 1 -^* by making 

Q h = h*Q€f2 1 (J 1 ir*) , n h = h*nef2 2 (J 1 ir*). 

Locally, the Hamiltonian section h is specified by h(t,q l ,pi) = (t,q l , —H,pi), where H is the 
local Hamiltonian function given by H = p^FC^Yv 1 — {FC~ V )*L. The local expressions are 

Q h = - Hdt , Q h = dq { A dp { + dH A dt . 

Of course J-C*@h = @£> an d J-C*Qh = 

The Hamiltonian problem consists in finding sections of f 1 , tp : R 5- J 1 ^*, characterized by 

^* i{X)Q h = , for every X € £( J 1 it*) . 

This condition leads to the Hamilton equations which, if ip(t) = (t,q i (t),pi(t)), in natural coor- 
dinates are 



dq i _ dH 
dt dpi 



dpi dH 
~dt ~ ~~dq J 



Assuming that these sections are integral curves of vector fields € X( J 1 ^*), the corresponding 
equations for these vector fields are 

i(x h )n h = o , i(x h )dt = \. 

As a final remark, it can be proved that solutions to the Lagrangian and Hamiltonian prob- 
lems are equivalent, in the sense that they are jF£-related; that is, 

y, = JF£ojV ; TJF£ o X c = X h o FC . (3) 



For regular, but not hyper-regular systems, the results are the same, but only locally on 
open neighbourhoods at every point, instead of J l ir*. 

A singular Lagrangian C is almost-regular if: V = TC^J 1 ^) is a closed submanifold of J l n* 
(let j: V J l n* be natural embedding), J Z C is a submersion onto its image, and for every 
y € J 1 ^, the fibres J-C~ 1 (J-C(y)) are connected submanifolds of J 1 tt. 

If C is an almost-regular Lagrangian, the submanifold V of J l ir* is a fibre bundle over E 
and M. In this case the /x-transverse submanifold j: V T*E is diffeomorphic to V. This 
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diffeomorphism is denoted by p: V >■ V, and is just the restriction of the projection u. to V. 

Then, taking the Hamiltonian section h = jo p _1 , we define the forms 

e° h = h*e ; n° h = h*n 

which verify that TCqQ^ = @c and FCqQ^ = &c (where TCo is the restriction map of TL 
onto V). Then, the Hamiltonian problem and the equations of motion are stated as in the 
hyper-regular case. Now, the existence of a solution to these equations is not assured, except 
perhaps on some submanifold of V, where the solution is not unique, in general. 

3 Unified formalism 

We define the extended jet-momentum bundle W and the restricted jet-momentum bundle W r 

W = JV x e T*E , W^AxgjV 

with natural coordinates (t,q l ,v l ,p,pi) and (t, q l , v l ,Pi), respectively. Natural submersions are 

pi : W JV , p 2 : W T*E , p E :W E , p R : W R (4) 

p\: W r , p r 2 : W r >JV , p r E :W r *E , p r R :W r >■ M. , 

with ir 1 o p 1 = t 1 o p o p 2 = p E . For y € J 1 -^, p G T*^, and [p] = p(p) € J 1 -^*, there is also the 
natural projection 

P w : W > W r 

(y,p) ^ (y. [p]) 

The bundle W is endowed with the following canonical structures: 

Definition 1 i. T/ie coupling 1-form in W is the p m -semibasic 1-form C G f2 1 (W) defined 
as follows: for every w = (j 1 ^), a) G W (i/ia£ is, aeT* («,)^ anc ^ ^ e T^W, i/ten 

c(y) = a(T U) (</,o /9R )y). 

,2. The canonical 1-form 6yv G J? 1 (W) is i/te p E -semibasic form defined by 0yy = P2@- 
The canonical 2-form is Oyy = -dByv = pffi G i? 2 (W). 

Being C a p E -semibasic form, there is C G C°°(yV) such that C = Cdt. Note also that 0>v is 
degenerate, its kernel being the /^-vertical vectors; then (W, fiyv) * s a presymplectic manifold. 

The local expressions for 0yy, Qyy, and C are 

@w = Pidq 1 + pdi , fiyv = — dpj A dcf — dp A dt , C = (p + PiV % )dt . 

Given a Lagrangian density C G i7 1 (J 1 7r), we denote C = p\C G i? 1 (W), and we can write 
C = Ldt, with L = p^L G C°°(yV). We define a Hamiltonian submanifold 

Wq = {w G W | £(u>) =C(w)}. 

So, Wo is the submanifold of W defined by the regular constraint function C — L = 0, which is 
globally defined in W using the dynamical data and the geometry. In local coordinates it is 

C-L=p + p i v i -L{t,q j ,v j ) =0 . 
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The natural embedding is jo : Wo W. We have the projections (submersions), see diagram (0): 

p\: W > JV , p° 2 : Wo T*E , p° : W *E, p^.Wo *- R 

which are the restrictions to Wo of the projections ([!]), and 

p° 2 =pop° 2 : Wo -JV. 

Local coordinates in Wo are (t,q l ,v l ,pi), and we have that 

Pi(t,q\v\pi) = {t,q l ,v l ) , 3o(t,q\v\pi) = (t, q\v\L — v l pi,pi) 
p%{t,q\v\pi) = (t,q\pi) , p%{t,q\v\p i ) = {t,q\L-v i pup i ). 

Proposition 1 Wq is a 1-codimensional p w -transverse submanifold ofW, diffeomorphic to W r . 



(Proof) For every (y,p) G Wo, we have L(y) = L(y,p) = C(y,p), and 

(a» w °Jo){y,p) = p w {y,p) = {v,Kp)) ■ 

First, p w o jo is injective: let (yi,pi), (?/2jP2) G Wo, then we have 

(Mw ° Jo)(yi,pi) = O w ° jo)(y2,Pz) =>• (m,p(pi)) = (y~2,p(p2)) =>yi = y2, m(pi) = m(p2) 

hence = £(2/2) = C(y\,pi) = C(y 2 , P2)- In a local chart, the third equality gives 

p(Pi) +Pi{Pi)v l {yi) =p{P2) +Pi(P2y(m) 

but/i(pi) = p(p 2 ) implies pi(pi) =pj([ Pl ]) = Pi ([p 2 ]) = p;(p 2 ); thenp(pi) =p(p 2 ), and pi = p 2 . 

Second, p w ojo is onto, then, if (y, [p]) G W r , there exists (y, q) G Jo(Wo) such that [q] = [p]. 
In fact, it suffices to take [q] such that, in a local chart of J 1 ^ x e T*E = W 

K(q) = Pi([p\) , p(q) = - p*([p]K(y) • 

Finally, since Wo is defined by the constraint function C — L and, as ker /iyy* = locally 
d ~ 

and — (C — L) = 1, then Wo is /x w -transversal. ■ 
op 

As a consequence of this result, the submanifold Wo induces a section h: W r *- W of the 

projection p w . Locally, h is specified by giving the local Hamiltonian function H = —L -\-piV % ; 
that is, h(t, q l ,v l ,pi) = (t,q l ,v l , —H,pi). In this sense, h is a Hamiltonian section of p w . 

So we have the following diagram 

(5) 




J L TT 
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Remark: Observe that, from the Hamiltonian |U w -section h: W r »- W in the extended 

unified formalism, we can recover the Hamiltonian ^-section h = j o p^ 1 : V >■ T*E in the 

standard Hamiltonian formalism assuming that C is almost-regular. In fact, given [p] G J 1 ^*, 
the section h maps every point (y, [p]) G (/>2) -1 ([p]) P^ [P2(h(y,[p]))]. Now, the crucial 

d 

point is the projectability of the local function H by p 2 . However, — - being a local basis 

ov l 

dL 

for ker p 2 *, H is /92-projectable if, and only if, pi = 7-^, and this condition is fulfilled when 

ov l 

[p]eP = ImFC C JV, which implies that p2[h((p r 2 )~ 1 ([p]))] G ^ = Im^£ C T*£. Then, 
the Hamiltonian section /i is defined as 

~h([p}) = (p 2 o hW2rHj(\p}))} = (jo A _1 )([P]) , for every [p] G P. 
So we have the diagram 




For (hyper) regular systems this diagram is the same with V = ImFC = J 1 ^* . 
Finally, we can define the forms 

e = f e w = P° 2 *e^^ 1 (m) , n = f n w = p° 2 *n£n 2 (w ) 

with local expressions 

9 = (L - piV l )dt + p i dq i , O = d(p i v i - L) A dt - dpi A dq* (6) 
and we have the presymplectic Hamiltonian systems (Wo,^o) an d (W r ,Q r ), with £l r = /i*$7o- 

4 The dynamical equations for sections 

Now we establish the dynamical problem for the system (Wo,^o) which, as a consequence of 
the diffeomorphism stated in Proposition [TJ is equivalent to making it for the system (W r ,O r ). 

The Lagrange-Hamiltonian problem associated with the system (Wo,^o) consists in finding 
sections of p®, tpo : M *• Wo, which are characterized by the condition 

r i(Y )Q = , for every Y G X(W ) . (7) 

This equation gives different kinds of information, depending on the type of the vector fields Yq 
involved. In particular, using /^-vertical vector fields, denoted by we have: 

Lemma 1 IJYqE is piL-semibasic. 

(Proof) A simple calculation in coordinates leads to this result. In fact, taking {^r} as a local 
basis for the /S^-vertical vector fields, and bearing in mind §§§ we obtain the /?^-semibasic forms 
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As an immediate consequence, when Yq G X V ^(Wo), condition (|7j) does not depend on the 
derivatives of ipQ: it is a pointwise (algebraic) condition. We can define the submanifold 

m = {(S,P) € Wo I i(V )(n ) { y, p) = 0, for every V G V (5ip) (p°)} 

where V^p®) denotes the /Og-vertical vectors. Wi is called the first constraint submanifold of the 
Hamiltonian pre-multisymplectic system (Wo, f2o), as every section tfto solution to (J7|) must take 
values in Wi . We denote by ji : Wi Wo the natural embedding. 

dL 

Locally, Wi is defined in Wo by the constraints pj = — -. Moreover: 

ov l 

Proposition 2 Wi is i/ie yrap/i of TL; that is, Wi = {{y,FC{y)) G W | y € J 1 ^}. 

{Proof) Consider y G J 1 ^, let R >■ be a representative of y, and p = FC(y). For every 

£/ G T^i^M, consider V = T^i (£)(/>(£/) and its canonical lifting V = T^ify\j (ft(U). From the 
definition of the extended Legendre map © we have {T yit 1 )* (F £{y)) = (Q c )y, then 

Furthermore, as p = J-£(y), we also have that 

i(V)[{Ty^y{TC{y))\ = i(T ff i (5) jV(^))[(T 5 7r 1 )*p] = i{{T ,^\{T ^^j 1 cft{U)))p 

= i(T^ {9) <f>(u))p = i(v)p. 

Therefore we obtain 

i(uw P ) = i(u){(j^r(e c )y) 

and bearing in mind the definition of the coupling form C, this condition becomes 

i(U)(C(y,p)) = i(U)[(j 1 cf>y@c)y}. 

Since it holds for every U G T^i^R, we conclude that C(y, p) = [(j 1 (ft)* © c]y > or equivalently, 
C(f7, p) = L(y, p), where we have made use of the fact that 0£ is the sum of the Lagrangian 
density C and a contact form i(V)dC (vanishing by pull-back of lifted sections). This is the 
condition defining Wo, and thus we have proved that (y,.F£(y)) G Wo, for every y G J 1 ^; that 

is, graph TL C Wo- Furthermore, graph TL and Wi are defined as subsets of Wo by the same 



local conditions: pi — — — = 0. So we conclude that graph T C = Wi- ■ 

As Wi is the graph of J-C, it is diffeomorphic to J 1 ^. Every section ipQ : R *- Wo is 

of the form ipQ = (iftc^n)^ with iftc = ° ifto'- K s- J 1 ^, and if V'o takes values in Wi 

then = T "£ o ^£ : R T*i?. In this way every constraint, differential equation, etc. in 

the unified formalism can be translated to the Lagrangian or the Hamiltonian formalisms by 
restriction to the first or the second factors of the product bundle. 

However, as was pointed out before, the geometric condition ([7]) in Wo, which can be solved 

only for sections ifto : R >■ Wi C Wo, is stronger than the Lagrangian condition tft* c i(Z)Qc = 0, 

(for every Z G ^(J 1 ^)) in J 1 ^, which can be translated to Wi by the natural diffeomorphism 
between them. The reason is that, as p\ is a submersion, and Wi is a p9-transversal submanifold 
of Wo (as a consequence of Proposition [2]), we have the splitting ^TW = TWi © Wl f^iPi), 
ji : Wi > Wo being the natural embedding. Therefore the additional information comes from 
the /^"Vertical vectors, and is just the holonomic condition. In fact: 
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Theorem 1 Let ipo : R >- Wo be a section fulfilling equation if) = (tpcipn) = (ipc^C > 

tp/l), where tpc = Pi ° V'o • Then: 



1. tpc is the canonical lift of the projected section (f) = p® E o ip : 
holonomic section). 



E (that is, tp£ is a 



2. The section ific = j 1 ^ is a solution to the Lagrangian problem, and the section p o ip-^ 
H o TC o ipc = o is a solution to the Hamiltonian problem. 



E such that j 1 ^ is a solution to the Lagrangian problem 



Conversely, for every section <p : 
(and hence J-C ° j l <fi is a solution to the Hamiltonian problem) we have that the section tpo 
(j'V, TL o j l cj)), is a solution to (?|). 

{Proof) 1. Taking |^-| as a local basis for the p\- vertical vector fields: 

d 



so that for a section i/jq we have 



and thus the holonomy condition appears naturally within the unified formalism. So we have 

(dq i dL \ 
t,q 1 ,— — ,— — I, since ipo takes values in Wi, and hence it is of the form ipo = 
dt ov 1 J 

(ffaFCofcj)), for <t>= {t,q l ) = p° E o^o- 
2. Consider the diagram 




T*E 



Since sections ipo '■ ^ *■ W) solution to take values in Wi , we can identify them with 

sections tp\ : R >■ Wi- These sections tpi verify, in particular, that ip* i(Yi)Oi = holds for 

every Y\ G 5E(Wi). Obviously ipo = Ji ° ipi- Moreover, as Wi is the graph of TC, denoting 

by p\ = Pi ° Ji '■ Wi >■ J 1 ^ the diffeomorphism which identifies Wi with J 1 -^, if we define 

ill = J i^o, we have that f2i = p];*J)£. In fact; as {p\)~ l {y) = (y,TC(y)), for every y G J 1 -^, 
then (pg o j x o (p^) -1 )^) = ^X(y) G T*£, and hence 

^ = o ji o (ply'yn = [((plr 1 )* o j! o p°*]o = [((pi)- 1 )* o j*]r2 = ((pir 1 )*^ . 
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Now, let X G X^ir). We have 



(p? o V> )* = (p° o Jx o Vi)* ;(X)^ 



(8) 



where Yq G 3£(Wo) is such that 1q = Ji*Yi. But as V'o i(^b)^o = 0, for every Yq G X(Wo), then 
we conclude that (j V)* z'(X)fi£ = 0, for every X G X(J 1 tt). 

Conversely, let j 1 ^: R >- J 1 ^ such that (j 1 ^)* «(X)f2c = 0, for every X G X(J 1 7r), and 

define V'o : R >- Wo as V'o = O'V) •^"-^ ° iV) (observe that Vo takes its values in Wi). Taking 

into account that, on the points of Wi, every Yo £ X(Wo) splits into Yo = Yq + Y 2 , with 
Y X G X(W ) tangent to Wi, and Y 2 G X V(p i ) (W ), we have that 

V> * i(Y )n = ^ i(Xo)^o + V'o ip^fio = 
since for Yq 1 the same reasoning as in (|8|) leads to 

VSi(Y 1 )o = (iV)*i(^o 1 )^ = o 

(where = (p})*^ 1 ), and for Y 2 , following the same reasoning as in 



a local calculus gives 



r i(Y 2 )n = (j^r 



fi(x) v, 



a 



dq l 
dx^ 



dt 



since is a holonomic section and Y 2 = fi — — . The result for the sections V'w = J~£ o j l cf) is 

OPi 

a direct consequence of the first equivalence relations ((3j) . ■ 
Remark: The results in this section can also be recovered in coordinates taking an arbitrary 

Oj r\ r\ rs 

local vector field Y = f— + f— +g l Tr-+ hi — G X(W ), then 

at dq l dv 1 dpi 



( dL 
i(Y )n = -f(p i dv i + V i d Pi -—Aq i 

-f l {^ + ^)+9 l {p, 



dL 



dv 
dL 

dv 1 



-6.V 



dt + hiiv'dt - dq i 



and, for a section V'o fulfilling ([ZD, 



= Vo i(Y )n 



dpi dL \ 



r -T- o2 - +g l ft-o-i 



(it dif / 



dv 1 J 



hi[v ~lt 



dt 



(9) 



reproduces the holonomy condition, the restricted Legendre map (that is, the definition of the 
momenta), and the Euler-Lagrange equations. The coefficient of the component / vanishes as a 
consequence of the last equations. 

Summarizing, the equation ([7]) gives different kinds of information, depending on the verti- 
cally of the vector fields Yq involved. In particular, we obtain equations of three different classes: 

dL 

1. Algebraic (not differential) equations, in coordinates pi = — t, which determine a subset 

Wi of Wo, where the sections solution must take their values. These can be called pri- 
mary Hamiltonian constraints, and in fact they generate, by projection, the primary 
constraints of the Hamiltonian formalism for singular Lagrangians, i.e., the image of the 
Legendre transformation, FC^tt) C J 1 -^*. 
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2. The holonomic differential equations, in coordinates v = , forcing the sections solution 

dt 

ipo to be lifting of 7r-sections. This property reflects the fact that the geometric condition 
in the unified formalism is stronger than the usual one in the Lagrangian formalism. 

3. The classical Euler-Lagrange equations, in coordinates 

d fdL\_ d 2 L d 2 q> | d 2 L dqi | d 2 L _ dL 
dt \dv l J dv^dv 1 dt 2 dqidv 1 dt dtdv 1 dq % 

dp ' d Ij 

which are obtained from — = — - using the previous equations. 

dt oq l 

5 The dynamical equations for vector fields 

Proposition 3 The problem of finding sections solutions to ffi is equivalent to finding the 
integral curves of a vector field Xq £ X(Wo), which is tangent to W\ and satisfies that 

i(x )n = o , i(x )dt = i. (ii) 

(Proof) In a natural chart in Wo, the local expression of a vector field Xq £ jE(Wo) is 

d ■ d ■ d 
^° ^ dt ^ dq l ^ dv l ^ % dpi ' 



Then, the second equation (|TTj) leads to / = 1, and the first gives 

coefficients in dpi : F 1 = v % (12) 

dL 

coefficients in dv l : pi = (13) 

av % 

dL 

coefficients in dq l : Hi = -^-r (14) 

dL f dL \ 
coefficients in dt : -F i — + G i [ Pi - ^ \ + Hy = . (15) 

dq l \ dv l J 

do dv l 

Now, if tp = (t, q t (t),v l (t),pi(t)) is an integral curve of Xq, we have that F % = — — , G % = ——, 

Lib lib 

dpi 

Hi = — — , and then (see equation 



• Equations (|12p are the holonomy condition. 

• The algebraic equations (|13p are the compatibility conditions defining Wi- 

• Using (|12|) and (|13|) . equations (|14|) are the Euler-Lagrange equations (jlOj) . 

• Taking into account (|12p and (|14|) . equation (|15|) holds identically. 

Observe that the condition that Xq (if it exists) must be tangent to Wi holds also identically 
from the above equations, since 

- X ( - — } - - G i - - v { + — (on W 

\ dv l ) dv l dvi dtdvi dq % dvi V dqi 1 

are the Euler-Lagrange equations again. Observe that, if L is a regular Lagrangian, these equa- 

dv i 

tions allow us to determine the functions G % = . If L is singular, then a constraint algorithm 

dt 
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must be used in order to obtain a final constraint submanifold Wf (if it exists) where consistent 
solutions exist, that is, Xq must be tangent to W/ (see [2] and Section [6] for details). ■ 

Now, the equivalence of the unified formalism with the Lagrangian and Hamiltonian for- 
malisms can be recovered as follows, where Xyvi(Wo) is the set of vector fields on Wo with 
support in Wi. 

Theorem 2 Let Xq be a vector field in Wo which is the solution to the equations Ml}) - Then 
the vector field Xc £ X(J 1 7r) ; defined by X^o p\ = Tp\oXo, is a holonomic vector field solution 
to the equations flj). 

Conversely, every holonomic vector field solution to the equations (QP can be recovered in this 
way from a vector field Xq £ Xyvi(Wo). 

{Proof) Let Xq be a vector field on Wo, which is a solution to (jlip . As sections ipQ : M. >■ Wo 

solution to the geometric equation ([7|) must take value in Wi, then Xq can be identified with 

a vector field X\: Wi >• TWi (i.e., Tji o X\ = Xq]^), and hence there exists X^: J 1 ^ 

>- T(J 1 7r) such that X\ = T(p\)~ 1 o Xc £ X(Wi). Therefore, as a consequence of the 

item 1 in Theorem [TJ for every section tpQ solution to (J7|), there exists X^ £ X(j 1 (j)(M.)) such 

that Tj^ o X® c = Xjc\ji^(v\, where j^: j 1 (/)(]R) 9- E is the natural embedding. So, Xc is 

^-transversal and holonomic. Then, bearing in mind that j\^Iq = p\*Qc> we have 

jt i(x )n = »(jfi)(j;no) = i(Xi)(pi*n £ ) = p\* i(x c )n c 

then i(Xc)Qc = because i(Xo)Qo = 0. A similar reasoning leads us to prove that, if i(Xo)dt = 
1, then i(X c )dt = 1. 

Conversely, given a holonomic vector field Xc, from i(X/;)f2£ = 0, and taking into account 
the above chain of equalities, we obtain that i(Ao)f2o £ [X(Wi)]° (the annihilator of X(Wi)). 
Moreover, Xc being holonomic, Xq is holonomic, and then the extra condition i(Yo) i(Xo)f2o = 
is also fulfilled for every Y £ X V(p?) (W ). Thus, remembering that jJTWo = TWi © Wl J*V(p?), 
we conclude that i(Xq)Qq = 0. To prove that if i(Xc)dt = 1, then i(Xo)dt = 1 is trivial. ■ 

Finally, the Hamiltonian formalism is recovered using the second equivalence relations ([3]). 
The proof for the almost-regular case follows in a straightforward way. 



6 An example: spatial semidiscretization of the nonlinear wave 
equation 

Consider the nonlinear wave equation given by 

d f do \ dg 

Utt = Tx\J^^ Ux) )-^ u) (16) 

where u : U C M? — > R, u(t, x) and a and g are smooth functions and we impose periodic 
boundary conditions u(t,x) = u(t,x + K), K > 0. Different choices of the functions a and g 
idealize one-dimensional models for fluids and materials. 

Equation (|16|) corresponds to the Euler-Lagrange equation derived extremizing the action 
functional 

(\ut -v(t,u x )-g(t,u)\ dtdx, 

where we will assume in the sequel the regularity condition — — 7 ^ 0. 

out 
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One basic idea towards a geometric discretization [SJ [TU\ [T3] of this type of equations is first 
to introduce an spatial truncation, that reduce the PDE (|16p to a system of ODEs preserving 
many of its geometrical properties. Hence, we replace the x-derivative in the Lagrangian by a 
simple difference (for simplicity, we will work with a uniform grid of N + 1 points, h = K/N) 
as follows: 

AT-l 

L(t,Ui, {ui) t ) = ^ 

i=0 



1 f (ui)t + (u i+ i)t 
2 



<T\t, 



Ui+1 



Hi 



9\t 



Ui+l + Ui 



In a more convenient notation, we are working with the Lagrangian function L : 
TR N+1 — ► M: 



N-l 



L(t,q\v*) = Y, 



i=0 



1 ^^-.(tyj-^i,^ 



where w J = and Q i+1 ' 2 = , i = . . . , N — 1, Q = q,v. Now, following the 

notation in previous sections, we find that 



9 = (L(t, q\v l ) - py) dt+p i dq i ,0 <i < N. 
Consider now a vector field 



satisfying the equations: 



It is easy to deduce that: 
f / = 1 

1 <9cr 

-no 



d d d d 

ot at/ 1 cw 1 opi 



1 



t,q 



0+1/2 



/i 9itx v 1 2 du 
1 f da i da 4 _j 



if 



N 



h \du x ' <9u 



1 /cfc 

2 1 



(*.9 i+V2 ) + £ (*.? i_1/2 ) 



/i 5u x 2du 
and the constraints defining VVi: 

1 



1 



Po = 2^° +1/2 > Pi = ^- 1/2 + pn = ^ 1/2 - 



Since Xq must be tangent to Wi then we obtain the additional conditions 



1< i <N- 1 



= X (po - ^ 0+1/2 ) 



= Xoipi-hv^ + v^ 2 )) 



,0+1/2 



G° + G 1 19a, o, 1 do , 



ft, V 



<9u. 



N-l/2\ 



4 ftdu/' j 2duV' 9 
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From these last equations we obtain G°, G 1 , ■ ■ • , G N 1 in terms of G N and the additional 
constraint 

£ <-!>'£ <«.-•) = o 

which determines the new constraint submanifold, W2. Again, the condition of tangency of Xq 
to W2 gives us a new constraint: 



JV-l 



i=0 



d 2 a , v t+l - v % d 2 a , , N 



du r dt 



0. 



determining the constraint submanifold, W3. From it, we obtain that 

JV-1 



8=0 



d 3 a , j\ v l+1 - v l d 3 a , A 



du T d 2 t 



+ 



Q3 



a 



du 3 



Qi+l _ Qi Q2 a 

h dui 



(t,vf) 



8 2 a 



which uniquely determines the remaining coefficient G form the regularity condition — — ^ 7^ 

out 



7 Conclusion and outlook 



Following the Skinner- Rusk model for autonomous mechanical systems, we have presented a gen- 
eralized framework for describing both Lagrangian and Hamiltonian time dependent mechanical 
systems. 

The key tool of this construction is the coupling form which is defined using the natural 
geometric structure of the manifold VV = J 1 ^ x g T*E, This function allows us to define in a 
natural way a submanifold Wq of W, which is diffeomorphic to W r = J 1 ^ x # J l n* , the true 
space of physical variables. Then, the compatibility of the dynamical equations stated in Wo 
gives a new submanifold Wi which is identified with the graph of the Legendre map TC, where 
all the characteristic features of the Lagrangian and Hamiltonian formalisms of time-dependent 
regular and singular non-autonomous systems are recovered. 

This unified formalism constitutes an alternative but equivalent approach to that given by 
Cantrijn et aim [2]. The essential difference is that, in this work, the dynamical equations are 
established directly in W = J 1 ^ x # T*E. These equations are compatible in a 1-codimensional 
submanifold of W, but the dynamical solution is undetermined, even in the regular case. In 
order to overcome this trouble, the authors are forced to introduce a new constraint, in such 
a way that the resulting submanifold is the graph of the Legendre map. As a consequence, 
they are unable to define intrinsically the submanifold of physical states Wo- In our model, the 
introduction of the coupling form gets round all the above problems. 

The Skinner-Rusk unified formalism which is developed here has been used to give a new 
geometric framework for time-dependent optimal control problems in PQ, where some interesting 
examples are analyzed. Following the above example the developed formalism could be applied 
to optimal control problems in partial differential equations where the spatial semidiscretization 
is used to solve. 
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